load datos_aeroCBA.dat

datos = datos_aeroCBA(:,1:2); 	% todas las filas de las columnas 1 y 2
filtro = ~isnan(datos(:,2));	% predicado que indica si existe valor en la i-esima fila de la columna 2

xanos = datos(filtro,1);		% filtra las x que cumplen filtro
xdatos = datos(filtro,2); 		% filtra las y que cumplen filtro

X = 1954:1:2013;
Y = interp1(xanos, xdatos, X,'spline');
plot(X,Y, xanos,xdatos,'r*');


